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Abstract 

This paper addresses the question of existence of (not necessarily 
self-similar) solutions to the 4- vortex problem that lead to total or par- 
tial collision. We begin by showing that energy considerations alone 
imply that, for the general A''-vortex problem, the virial being zero is 
a necessary condition for a solution to both evolve towards total col- 
lision and satisfy certain regularity condition. For evolutions assumed 
to be bounded, a classification for asymptotic partial collision config- 
urations is offered. This classification depends on inertia and vorticity 
considerations. For non-necessarily bounded evolutions, we discuss the 
relationship between partial and (non-regular) total collisions when the 
virial is not zero and a generic condition on the vorticities holds. Fi- 
nally, we give a canonical transformation that, for a class of 4- vortex 
systems satisfying a restriction on the vorticities, allows to formally 
apply the averaging principle in order reduce the dynamics when two 
of the vorticities are near a binary collision. The reduced system is a 
one-degree of freedom hamiltonian system. 

1 Introduction 

The study of the dynamics of point-vortices on two dimensional manifolds, 
for an ideal incompressible fluid, has become a very active field of research 
in recent years. For background on the subject we refer the reader to the 
comprehensive book of P. Newton (2000) j9j. The case of vortices on the 
plane has a long history going back to W. Grobli's 19th century work (cf. 
0]), who studied the analytical solutions of the three vortex problem for 
certain special cases. Further developments on the three-vortex problem 
were carried out in ^ El El E] • In |H| (see also jlj ) we prove that total 
collision for three point vortices on the plane is only possible when the virial 
(defined below) is zero, in which case the vortices evolve in a self-similar 
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fashion. In op. cit. we also prove that partial (i.e. binary) collision can not 
happen with only three vortices on the plane. 

In this paper we are interested in the problem of existence of total or 
partial collapse for a system of four point-vortices on the plane. The motion 
of four vortices in an arbitrary configuration, but all with equal vorticity, has 
been studied in |^ (see also both analytically and numerically, where 
it is shown that the four-vortex problem is not integrable. In that paper 
a canonical transformation reducing the system from equal vortices to a 
system with N — 2 degrees of freedom is found. A similar transformation 
is shown to exist in [7] for the dynamics of two vortex dipoles where the 
total vorticity of each dipole is zero and an analytical study of collision 
trajectories and a numerical study of the scattering in both the regular and 
the chaotic regime is offered. 

To our knowledge, all collisions solutions to the A^-vortex problem stud- 
ied in the literature involve self-similar evolutions. It is known (^J) that 
these can only occur when the virial is zero. So a first natural question 
is whether a null virial is a necessary condition for total collision. We an- 
swer this question affirmatively for a class of collisions that are regular in 
a precise sense. To illustrate this result, we apply it to the problem of four 
vortices in parallelogram configuration, first studied by Novikov and Sedov 
in op. cit. After defining carefully the different senses in which a system 
can develop partial collisions, we determine necessary conditions for partial 
collapse and necessary asymptotic configurations under the assumption of a 
bounded evolution. Then, for non-necessarily bounded evolutions we study 
the relationship between partial and (non-regular) total collision, classify- 
ing this relationship in terms of the sign of the vorticity. In this context a 
non-zero virial and a generic condition on the vorticities is assumed. 

Some of the results described in the previous paragraph motivate the 
study of the evolution of a 4-point vortex system in a configuration where 
two of the vortices (say Ti and r2) are close to binary collapse. This is a 
complicated problem but we believe that averaging methods can give useful 
information. A first step in this direction is achieved in the last part of 
this paper, where we describe a canonical transformation that, in the class 
of systems satisfying the condition Fi -|- r2 = = allows to formally 
apply the averaging principle in order to reduce the 4-vortex system to a 
one-degree of freedom hamiltonian system. 
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1.1 Preliminaries 



The system of N point-vortices on the plane can be described as fohows. 
Let = Xq + iT/o- G C be the coordinates of the a-th point vortex with 
vorticity T^, a = 1, . . . , N. Thus the configuration space is C" = M^". The 
equations of motion are given by 



N 
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This system is Hamiltonian, with the energy (Hamiltonian) function given 
by 

H = -^J2^-'^P^^^o.p, (2) 



2tt 

a<f3 



where l^p := \za — zp\, and symplectic form given by 

n 

r2(z, -w) = - Im ^ TaZaWa , (3) 
a=l 

where the G M, Tq 7^ 0. Hence H is a conserved quantity, but also some 
other conserved quantities arise. From it follows that 

Z = ToZq, 

a 

and 

I = ^'^pI: Pa = \Za\ 
a 

are conserved as well. Here F = Fq. 

If F 7^ we call Z/T the barycenter of the system. A calculation shows 
that 

M:=^F„F^/2^ = F/-|Z|2. (4) 

a<l3 

Thus, conservation of Z and / imply conservation of M. Finally, from 
it also follows that the angular momentum or virial 

y = ^ F« Zq, X io, = ^ TaTi3 

a a</3 

is an invariant of the system, i.e., a constant depending only on the vortici- 
ties. Without loss of generality we will assume that Z = when F 7^ 0, and 
hence M = TI. 
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Notice that all initial configurations that end in total collision must sat- 
isfy M = I = (even when F = 0). 

We say that the dynamical evolution of N vortices on the plane is a self- 
similar trajectory if it can be expressed as Zi{t) = C(t)zi{0), i = 1, . . . N, 
for some ( : M — > C (hence Z = 0). 

1.2 Total and partial collapse: definitions and generalities 

In order to be able to state the results in later sections, we must be precise 
on what we mean by total or partial collision in the A^-vortex problem. (We 
remark that the following definitions are independent of the specific dynamic 
structure of point-vortex systems inherited from Euler's equation.) 

Definition 1. We say that a system of N point- vortices (with positions 
Zi ^ C, i = 1,...,N) evolves towards an n-collision (binary, ternary, 
etc. -collision if n = 2, 3, etc.), or simply a collision, in time < oo if there 
are indices / = {ii, . . . , n < N, such that for every i,j £ I we have that 
\zi—Zj\ ast ^ T^,. We say that the system evolves towards a sequential 
n-collision (binary, ternary, etc.), or simply a sequential collision, if there 
is a monotonically increasing sequence of times {tk} converging to T.^ < cxd 
such that \zi{tk) — Zj{tk)\ — > as A; — > oo for all i,j £ I. 

It has been shown in 8^ that in the three-vortex problem there are no 
binary collisions. However, it makes sense to consider the existence of partial 
collisions in the A^-vortex problem. To this end we state with the following 

Definition 2. If in the previous definition n = N we say that a system 
evolves towards a total collision. Assuming that the vortices labelled 
ii,...,in, n < N, evolve towards a total collision, we call it a proper 
n-collision (proper- binary, ternary, etc.), or simply a proper partial col- 
lision, if there are no other vortices labelled in+i, • • • , in+k^ k < N — n, such 
that the whole set of vortices ii, . . . , in+k evolve towards a total collision. A 
proper sequential n-collision is defined analogously. 

It also makes sense to talk about partial collisions in the context of 
similarity classes of A^ point-vortices on the plane. So, for example, it makes 
sense to speak of a binary collision within a process of total collision or blow- 
up. To define these types of collisions we first introduce a notation for the 
rescaling of the (squared) distance between the point-vortices: 

If {zi{t), . . . ,ZN{t)) is a trajectory in describing the evolution of a 
system of A^ point- vortices, let bij := \zi — Zj\'^ and let the parameters {p, Aj} 
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be defined by the conditions bij = pPij, p > 0, Yli<j (^ij — 1- With this 
notation we state the following 

Definition 3. We say that a system of N point-vortices evolves towards a 
relative n-collision (relative- binary, ternary, etc.), or simply a relative 
partial collision, in time T,,, < oo if there are indices / = 
n < N, such that for every i,j £ I we have that f3ij — > as t — > T*. A 
relative sequential n-collision is defined analogously. 

The term proper relative (sequential) n-collision is defined in an 
obvious way. 

Remark. It is clear that {proper n-collisions} C {proper relative n-collisions} C 
{relative n-collisions} C {n-collisions}. 

In this paper we consider the question of the existence of evolutions of 
the four vortex problem that lead to total collision when the virial is not 
zero. From the energy relation Q it follows that if total collision is to occur 
then the vorticities can not all have the same sign. The two possible cases 
are: a) three vorticities with the same sign and one with opposite and b) 
two positive and two negative vorticities. 

Let bij := \zi — Zj\'^. Then a quadrilateral configuration is represented 
by a point (612, 613, 614, 623, 624, 634) € M^, subject to the constraint that the 
Cayley-Menger determinant det(Vij) vanishes, where {Vij) is the 5x5 matrix 
given by 



where 5ij is the Kronecker delta, representing 288 times the square of the 
volume of the tetrahedron with sides of length \zi — Zj\ (cf. jl4l p. 256]). 

Through out this paper we use the rescaling coordinates {p, (3ij} defined 
by hj = pPij, subject to the constraint Pij = 1- Notice that to each 
point (Pij) £ satisfying both Pij = 1 and det(V^j) = we can asso- 
ciate a point in the space of quadrilateral shapes, understood as the set of 
equivalence classes in the set of all quadrilaterals, with equivalence given by 
similarity, and that this correspondence is a bijection. 

Exponentiating both sides of @ we arrive at 




f3ri^2 ffi^s ^^2^3 /prir4 ar2r4 ai 

pV ~ ^12 Pi3 P23 ^14 ^24 Pi 



(5) 
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2 Dynamics and total collapse for point- vortices 

2.1 Dynamics in the space of squared distances 

The discussion of the cohapse of N point-vortices on the plane is best dealt 
with by describing the dynamics in the space of squared-distances between 
the vortices. This space is given by the set Q of points (612, • • • , 6(7v-i)Af) G 
M^s) such that bij := \zi — Zj\^ , where (zi, . . . , zjy) G C^. Thus, Q is the 
set of vectors in MV2j with non-negative entries such that, for every triad 
{i,j,k) of vortex labels, {bij,bjk,bf:i) satisfies the triangle inequality 

bfj + bjk + bli < 2{bijbjk + hjbki + bjkbki) ■ 

Let C be the collision subset of Q. That is to say, C is the set of points 
{bi2, . . . , ^(Ar-i)Ar) £ Q such that bij = for at least one pair of vortex-labels 

The dynamics in Q is given by the following 

Proposition 4. For a system of point-vortices {Ti, . . . ,r7v} the equations 
of movement are given by 



2^Er.^./=.(^-^) , ^<J■ (6) 



Hk ^kj 



Here bij = \zi — Zj\'^ , Zi ^ C denotes de position of the i-th vortex and An^j 
is the oriented area of the triangle A{zi, z^, Zj). 

For the proof of this proposition we need the following 

Lemma 5. The identity 

Im ( — + —] =Im{ziZ2) 



_Z2 Zi/ ■ ■ \\Zi\'^ \Z2^'^ 

holds for all zi, Z2 £ C \ {0}. 

Proof. Express the Zi in polar coordinates. □ 
Proof of proposition^ Let Wij = Zi — Zj for i ^ j. Then 
; d , _ , . _ 

^ij = ^ [WijWij) = WijWij + WijWij . 



6 



From dU, 



Hence, 



Wi 



2tt 



Wik Wij 



k¥=i,j 



Wjk Wji 



2tt 



Wk wkj I 



and 



it. ^^^^ ""'^^ 



Summing both equations, while taking into account that 

Wjj _^ Wij_ _ 2 + !^ + ^ 

Wik Wkj Wik Wkj ' 



we obtain 



6,- = ^Er,Imf!^ + ^ 



k^ij 



Applying lemmajSland noting that Ai^j = Im {wi^Wkj) /2 we get the desired 
result. □ 



2.2 A criterion for total collapse of N point-vortices 

We have seen that a condition for total collapse in the point-vortex prob- 
lem is that M = 0. In this subsection we will show that, under the assump- 
tion that the evolution of the system satisfies certain regularity condition, 
another necessary condition for total collapse is that V = Q. As we will see, 
the proof of this fact is based solely on energy considerations. 

In what follows we will use the following standard notation. If w G 
M" \ {0} then v := v/\\v\\. If f{t) is differentiable then / := df/dt. 

Lemma 6. Let 7 : (a, b) — > he a regular curve not intersecting the 
origin, b < 00, such that limt^f,7(t) = 0. Then the following are equivalent: 

i. The limit of ^{t) as t ^ b exists. 

a. The limit of ^{t) as t ^ b exists. 
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Moreover, if these limits exist then one is the negative of the other. 
Proof. Suppose that i. holds. By L'Hopital's rule, 

lim7(t) = lim -P^^ = lim — . 

'-'mhm '-'lit) -at) 

The last equality follows from the identity ||7||J^||7|| =7'7, easily obtained 
by time-differentiating ||7|p. Since the limit in the left-hand side exists then 

1 



lim llT'ft)!! = lim ■ 



Hence limt^b^{t) ■ 'y{t) = ±1 and it follows that ii. holds. 

Suppose now that ii. holds. Let t : (a, c) (a, b) be a reparametrization 
such that 7(t) = 7(t(r)) is a curve parametrized by arc length. If c = co 
then the fact that limT-^oo 7'(''") exists, with ' = d/dr, implies that ||7(t)|| is 
unbounded as r — > oo, contradicting the hypothesis that limT-^oo 7('7") = 0. 
Thus c < oo and there exists curve parametrized by arc length 7 : (a, c+e) — > 
M", e > 0, such that 7(0) = and 7 = 7|(a,c). Then j{t) = -^'{c){c - 
t) + /o(t), where p{t) = o(c — r). Hence 

^ _^'(c) + gli 
lim7(t) = lim 7(r) = lim = —7 '(c) , 

t^b'^^ r-^c- ^ ' r^c- II _ ^ /(c) + £M II 

showing that i. holds. Finally, it is clear that 7 '(c) = \mit_,i,"f[t). Hence 
the limit in i. and minus the limit in ii. are equal. □ 

Definition 7. Let 7 : (a, h) — > M" be a regular curve not intersecting the 
origin, h < 00. We say that 7(t) approaches the origin in a regular 
way as t ^ 6 if limt^hj{t) = and 7 satisfies any (and hence both) of the 
properties i. and ii. in lemma El 

Remark. An analogous lemma and definition can be stated for a regular 
curve 7 : (a, b) — > R" when ^{t) — > as t — > a. It is easy to see that the 
only difference in this case will be that the limits in i. and ii. of lemma ^ 
will be the same. 

We now state the main result of this subsection: 

Propositions. Leth{t) = (^i2(i)> • • • > ^(Af-i)Af(t)) be the curve in'M.^'^' that 
is the solution to the equations of motion © for a given initial condition 
b(0) G Q\C. Suppose that there is a time < 00 such that the curve h{t) 
approaches the origin in a regular way as t ^T^,. Then V = Q. 
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Proof. The set Q \ C is contained in Q := the set of vectors in W^^) with 
positive entries. From @ we have that, on Q, the surfaces of constant 
energy are given by £h = F^^{h), where h = exp(— 47r//) and 

F:Q-^R: F(6i2, . . . , \n-1)n) ■= J] ^f' ■ 

i<j 

Let b = (612, . . . , ^{Ar_i)Ar) £Q, 5h = ((56/2, • • • , Sb(^N_i)isi) G TbQ. Since 

dF, . 5h = F(b) ( r,r2^ + • • • + r(,.,)r.^^(^) (7) 

then F is a submersion. Hence Eh is a smooth manifold whose tangent space 
at b is given by kerdFb. From ((7|) it is clear that, as subspaces of M^z)^ 
kerdFb = kerd-Fxb for all b G Q, A > 0. Therefore, 

where w := limt^T* b(t) is guaranteed to exist by the assumption that h(t) 
approaches the origin in a regular way. 

On the other hand, from lemma IHl and the remark following definition 

|71 tii = ±limt_+T, b{t). Since b{t) G Tf,(i-^£h for all t, we have that w G T^Sh- 
Hence, 

= dF^-w = F{b){TiT2 + ■■■ + r(^N~i)^N) ■ 
Thus V = 0, as claimed. □ 

2.3 Application: total collision for a vortex parallelogram 

In this section we consider a system of four point- vortices in parallelogram 
configuration with the property that the vorticities at opposite vertices are 
equal. (This problem is discussed in assuming from the start a self- 
similar collision.) We will use proposition |H1 to show that the system admits 
total collapse only if V = 0. 

Consider then a four point- vortex with Fi = Fs, r2 = F4 and an initial 
condition satisfying 612 = 634,613 = 624. Then, fromEl it follows that 

. 1 \ p ) , 

27r \b,k bjj ^ " ' 

with k, I) any permutation of (1, 2, 3, 4), and hence 612—634 = 613 — 624 = 
0. This shows that the property of being in a parallelogram configuration is 
preserved through out the evolution of the system. 
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Now let us consider evolutions leading to total collapse. Then we have 
to assume that M = and rir2 < 0. The conservation of energy H and 
inertia M, together with the parallelogram law 

bi2 + &34 = 2(6i3 + hi) 

imply that 

iriFsl (6i3 + 634) = r? 613 + Tl 624 , bll b^l = h (612 buf^'^'\ . (8) 

Noticing that the first equation in Q implies that 613 = |r2/ri|624, it follows 
that system Q is equivalent to 

h{P)z = h{xy)\ {l + p)z = 2{x + y), (9) 

where {x,y,z) = (612,614,624), P = 1^2/^11, 6 = 2|rir2|/(r2 + F^) and 
fi{x) = \ Eliminating z from Q we obtain 

f2{f3;h)ix + y) = {xy)' , 

where f2{x;h) = 2/i(x)/(/i(l + x)). Thus, introducing the projective co- 
ordinate p = y/x one gets the following reparametrization of the evolution 
curve of {x{t),y{t), z{t)), 

-i^' -(c^F' <^»' 

where A = h/fi{P), 7 = f2{(3;h)~^ and a = 1/(1 — 2(5). An analogous curve 
is obtained using the projective coordinate q = x/y as parameter. Notice 
that 1 - 25 = (iTil + iraDVdFip + jFap) > 1. Hence < a < 1. 

We are interested in the case when {x,y,z) — > as p ^ (or q 0). 
This implies that < 6 < 1, which we now assume in what follows. For 
p ^ 0, {x{p),y{p), z{p))^ = c{p) where 

c{p) = . 

Dividing c{p) by p"^ (by ^"(i-'^)) if < 5 < 1/2 (if 1/2 < 5 < 1) we obtain 
that limp_»o c(p) equals (1,0,0), (1,1,0) or (0,1,0) depending on whether 
< 6 < 1/2, 6 = 1/2 or 1/2 < 6 < 1, respectively. Hence {x{p),y{p), z{p)) 
is a curve that approaches the origin in a regular way as p — > 0. It is 
clear that an analogous conclusion is reached if we work with the parameter 
q = 1/p instead. Therefore, by proposition |H1 the four- vortex problem in 
parallelogram configuration under consideration admits total collapse only 
if F = 0. 
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3 Necessary properties of collision evolutions in 
bounded dynamics 



In this section we establish some necessary properties that any candidate 
dynamic evolution, leading to binary or ternary collisions in the 4-vortex 
problem, must satisfy. The properties to be discussed refer to the distances 
between the various elements of the limit configuration (as the time tends to 
the collision time) and necessary relations to be satisfied by the vorticities. 
These properties are obtained only from energy and inertia considerations. 

The energy equation @ alone restricts the types of collisions that can 
occur if we assume that the system remains bounded for all time. Indeed, 

Proposition 9. // a system of four point vortices remains in a bounded do- 
main while it evolves towards a collision then it is either a ternary collision, 
a double binary collision or a total collision. 

Proof. The energy equation El implies that if the system is bounded and 
develops a collision then not all the vorticities have the same sign. W.l.o.g. 
we can assume that either a) ri,r2 > and r3,r4 < or b) ri,r2,r3 > 
and r4 < 0. 

Case (a). Equation becomes 



If one of the factors in the l.h.s. of tends to zero (two vortices of the 
same sign collide) then at least one of the terms in the r.h.s. tends to zero 
(two vortices of opposite sign collide). This leads to a ternary collision. 
Similarly, collision of two vortices of opposite sign implies collision of two 
vortices of the same sign and again we have a ternary collision. Notice that, 
if M = 0, the ternary collision can be part of a total collision. 

Case (b). Equation © becomes 



Suppose that two vortices of the same sign collide, say bi2 — > 0. Then at 
least one of the factors in the r.h.s. of (|12|) must tend to zero. There are 
two distinct possibilities, represented by 614 — > and by 634 — > 0. (The case 
624 — > represents the same case as 614 0.) The first possibility leads 
to the ternary collision (1,2,4) and the second leads to the double binary 






(12) 
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collision (1,2), (3,4). If we start by assuming that two vortices of opposite 
sign collide then at least one of the factors in the l.h.s. of H12() must tend to 
zero and the same argument applies. □ 

Remark. In the context of the proof of the previous proposition, it is easy to 
check that in case (a) (i.e. Til, r2 > and T^, T4 < 0) there are four ways in 
which the vortices can conglomerate in a candidate proper ternary collision, 
namely, one for each vortex excluded from the collision of the other three. 
And in case (b) (i.e. ri,r2,r3 > and r4 < 0) there are three possibilities 
for a proper ternary collision, namely one for each positive vortex excluded, 
and there are three possibilities for a double binary collision, namely one for 
each of the positive vortices colliding with the negative one while the other 
two positive vortices collide. 

3.1 Conditions and configurations for candidate partial col- 
lisions 

Let us consider what are the implications of conservation of M for any 
candidate evolution leading to a (not necessarily proper) triple collision or 
double binary collision. All vortices in the system are assumed to have 
non-zero vorticity. 

Ternary collision. Suppose that there is an evolution leading to the col- 
lapse of vortices (1, 2, 3). Then 614 — 624 — *■ and 614 — 634 — > 0. Therefore 
Q implies that, in the limit as time approaches the collision time, 

M = (Ti + T2 + T3)T4(f (13) 

where = 614 = 11^123 "-^Hp. Here 2^23, z* are the positions of the colliding 
point of vortices (1, 2, 3) and of the limit point of vortex 4. 

Double binary collision. Suppose that there is an evolution leading to 
the simultaneous collapse of vortices (1,2) and (3, 4). Then Q implies that, 
in the limit as time approaches the collision-time, 

M = (ri + r2)(r3 + r4)d2, (i4) 

where (P = 613 = Hzjj'g, Here the positions of the colliding 

points of vortices (1,2) and (3,4). 

Prom the discussion of the cases leading to equations H13() and ((Tl]) we 
arrive at the following collision conditions depending on whether M is equal 
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to zero or not, on the vorticities and on the limit configuration at colhsion 
time: 

Proposition 10. Consider a system of four point-vortices each with non- 
zero vorticity. Use the notation Fjj = Ti + Tj and Tij/. = Fj + Fj + rfc. Then 
necessary conditions for ternary or double binary collapse are given by the 
following table: 



Collapse type / 
vortices colliding 


Condition on M 


Vorticities 


Limit configura- 
tion parameter 


Ternary 


ifM = 


^ijk = 


d arbitrary 




MFyfeF, >0 




Double binary 
{i,j) and {k, I) 


ifM = 


TijFfc; = 


d arbitrary 


ifM^O 


MF^^ Fm >0 





Here d is the distance between the limit cluster points, as time approaches 
the collision time. 



4 On the relationship between partial and total 
collapse 

In this section we drop the assumption of bounded dynamics. In the propo- 
sition and corollary that follows we consider the setting of a system o four 
point vortices, all with non-zero vorticities and such that two of them, say 
Fi and F2, are of opposite sign. 

Proposition 11. Suppose that {Ti+r2){Ti+r2+T3)-{ri+r2+T4) / and 
AI = 0. Let {l3ij} be the normalized square distances between the vortices, 
as defined above. Suppose that, as t ^ Tq < 00, P12 0. Then there is a 
T < Tq and there are constants < ci < C2 such that if T < t < Tq then 
ci < (3^12^^^^ / < C2, where V is the virial and p is the quadrilateral scale 
factor defined above. 
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Proof. Noticing that P12 — > implies that P23 — P13 and ^24, — Pi4 
we obtain, from ©, that 



h(3 



12 /ori2r3 /ori2r4 



where := Fi + The constraint ^ = 1 imphes that 

2/?13 + 2/3l4 + /?23 = 1 + 9^2 (t). (15) 
The condition M = imphes 

ri2r3/3i3 + ri2r4/3i4 + TsVipsA = ^sit) . (i6) 

Here the functions ipi{t), i = 1,2, 3, tend to zero as t ^ Tq. Since the triple 
(/5i3, /3i4, /334) represent the squares of the sides of a triangle then it lies in 
the cone 

Ph + Pu + PL < 2 iPisPu + P13P34 + PiM . (17) 

The intersection C of the line given by the limits, a t ^ Tq, of (|15j) and (|16)) 
with the cone (|T7|) . if not empty, is a closed segment and thus a compact set. 
(If >C = then P12 can not tend to zero and the proposition is trivially true.) 
Identify {P13, Pu, Psi) with point {x,y,z) G M'^. The cone (fT7|) is tangent to 
the cartesian planes at the lines {x = y,z = 0), {x = z,y = 0, {y = z,x = 0). 
The hypothesis (Ti + r2)(ri + r2 + Tg) • (Fi + r2 + r4) / implies that 
these lines do not intersect ()16p in the limit as t ^ Tq. Therefore xyz^O 
for all (x, y, z) £ C. Let 

(ci, C2) = (min, max) x^^^ra yViiTi ^rgr* _ 

{x,y,z)eC 

Since the distance between the point (x, y, z) and the segment C tends to 
zero as t — > To then, for any 5 > 0, there exists T < Tq such that if 
T <t <To then ci - 5 < x'^i^rs yri2r4 2r3r4 < C2 + 5. Choosing (ci, C2) = 
(ci — (5, C2 + 5)//i, with (5 = ci/2, we get the desired result. □ 

Corollary 12. Suppose that vortices 1 and 2 have vorticities of opposite 
sign, that (Ti + r2)(ri + r2 + r3)(ri + r2 + r4) / O and thatM = 0. Then, 
if V > (resp. V < 0), P12 — > only if p ^ 0, (resp. p — > ooj. Thus, 
a relative binary or ternary collision involving vortices 1 and 2 can occur 
only as part of a process of total collision ifV>0 and as part of a blow-up 
process ifV<0. 
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Proposition 13. Suppose that a system of four point vortices, all with non- 
zero vorticity, evolves (maybe sequentially) towards total collision. Then the 
system evolves towards a binary or ternary sequential collision. 

Proof. If the system evolves towards total collision then there is a sequence 
of times {tn}, tn ^ T^, < oo, such that {p{tn)} is a sequence converging to 
zero. It follows from Q that if y > (resp. V < then there must be 
at least one pair of labels with TiVj < (resp. TiVj > 0) such that 
{f^ijitn^)} is a sequence converging to zero, with {tn^} a subsequence of 
{tn}. This implies that vortices i and j evolve weakly towards collision. □ 



5 Dynamics of a four- vortex system with a binary 

We have seen that, when the virial is positive, total collapse has to be 
accompanied by a binary or ternary collapse. In this section we consider 
the dynamics of a four-point vortex system in a state very close to having 
a binary collision. The two vortices that are near collision can be called a 
vortex binary. 

Consider then a system of four point vortices with /3i4 << 1. Let its 
configuration be given by {zi, Z2, zs, Z4) G C^, so that \zi — Z4I is very small 
relative to (some of) the other distances. Introduce the linear change of 
coordinates 



Z2 = Z2, Zs = Zs 



Z4. 
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1 14 



with := Fi + g, C G C. Thus q is the vector going from zi to Z4 
and C is the position vector of the center of vorticity of the pair (1,4). A 
calculation shows that the vector field in the (g, (, Z2, Z3) coordinates is given 

by 



ri4-|^ + r2 {rI - R2) + Ts [rI - Rf) 



C = [r2 {TiRl + URl) + Fs {TiRl + URl)] 



Z2 



Z3 



1 

i 

2^ 



Ts 

r2 



Z2 - Z3 

\z2 - zal^ 

Z3 - Z2 
\Z3 - 



F4i?2 — Fii?2 

— TiRi 
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where 



= Fi + r4 and 



^ ~ + m^l I'' ^1 ri4 

j = 2,3. 

These can be regarded as the dynamic equations of two coupled systems, 
{q\ and {(", Z2, -23}. Notice that, while at (? = the dynamic equation for {q} 
becomes singular, the dynamic equations for {C, ^^25 -Zs} become the equations 
of a standard three-vortex system with vorticities (Fi + r4,r2,r3). 

Let us now consider closely the system {q}. With Wj = Q — Zj, j = 2, 3, 
the terms — Rj) appearing in the expression for q can be written as 



R]-R'' = - r ,0,^'' r ,o 1 + -, o, (18) 



jl d4 _ I -'I i 14 

Fi |2| 

where Wj ■ q := Re{wjq). Thus, 



1 = ^ [[j^ + /l(^^2,1f3,g)j g + f2{w2,W3,q)w2 + h{w2,W^,q)w^ 

(19) 

with /l(w2,'W3,0) = T2/\W2\'^ + r3/|'U;3p, /2(l«2,U^3,0) = /3 (w^2 , "^3 , 0) = 0. 

We can regard (|19|) perturbation of 

q = -^{^^ + h{w2,W2„q)^ q (20) 

whose phase portrait consists of concentric circles around the origin. 

We now study the change in the phase portrait due to the perturbation 
term i(/2'W2 + /3if3)/(27r). To this effect we need to understand some general 
aspects of the relationship between the dynamics of {q} and {(", Z2, z^}. 

The first observation is that the total energy h is approximately the prod- 
uct^ of the energies of {q} and {C, Z2, 2:3}. More precisely, given (zi, . . . , Z3) G 
C", let /i(zi, . . . , z^; Ti, . . . , r„), M(zi, . . . , z^; Ti, . . . , r„) denote the asso- 
ciated energy and the value of Ylii<.j '^i^j\zi — Zj\'^ , for the given vorticities 
(Fi, . . . , r„). It is easy to see that if \q\ « q (i.e. the binary is very small) 
then 

^(2:i,2;2,2;3,^4;ri,r2,r3,r4) h{Q, z2, z^-Tii,T2,T^)\q\'^^^^^ . 



^We talk about the product and not the sum because h is obtained from the exponen- 
tiation of the original energy introduced in 
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The second observation is that, if the vorticities satisfy certain inequahty, 
then for any given energy, if the binary is small enough, the dynamics of 
{C, Z2-, Z'i} is much slower than the dynamics of {q}. More precisely, 

Proposition 14. If V > 2|rir4| then given /iq > and e > there ex- 
ists 6 > such that if \q\ < 6 then, for any (C) -22)^3) £ such that 
h{C,Z2,Z3;Ti,T2,T3) = ho/\q\^^'^^ and M(C, ^2, ^3; Tm, Ts, Ts) = 0, we 
have that 

\q\ 

Sketch of proof Let r(t) = {bi{t),b2{t),h(t)), hi = |z2-Z3p, ^'2 = |C--Z3p, 
63 = IC — Z2p- We know that, denoting 7(t) = 7_R(t) + ^B{t), where 
7_R is the radial component of the velocity and 75 is the component par- 
allel to the plane B : 61 + 62 + 63 = 1, 7_B(t) = fo/||r(t)|| for some 
constant vector vq parallel to B. Also, if Pminih) denotes the minimum 
value of ||7(t)||, with j{t) constrained to have energy h and M = 0, then 

Pmm(/i(C,^2,^3;ri4,r2,r3)) ~ /i(c,^2,^3;ri4,r2,r3)V^ = ho\q\-^\^^^^\/^ . 

On the other hand, |g| ~ l/\q\- Therefore |7_b|/|9| ~ |g|i^2|rir4|/v 

result follows. □ 

Remark. We have seen above that it can happen that, besides /3i4, /3j4 tends 
to zero as t — > io for j = 2 or 3 (but not both). If such is the case then 
we have a ternary collision and thus \q\ and \wj\ may be comparable. One 
can study the case when this is not the case, so that r4 participates only 
in a binary collision. More precisely, one can assume that << 1, 

j = 2, 3. In this case the denominator in both terms of the r.h.s. of (|18|) can 
be approximated by \wj\'^ and the perturbation term can be approximated 



as 



^(/2u;2 + /3u;3)-^ 



^I24\q\ - 2w2 ■ q , ri34 - 2'u;3 • g 

, [4 W2 + . 71 W3 



with Tij4 := rj(r4 — ri)/(ri + r4) and with 
rij4kP - 2ii;j • q 



Wj«l, j = 2,3. (21) 



5.1 A canonical transformation allowing averaging 

We have seen above that the equations of motion for the four vortex problem 
with a binary can be thought of two coupled systems, one for the binary 
and one for a three-vortex system, with an added perturbation that couples 
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the system. In this subsection we wih discuss a canonical transformation 
that ahows to formally apply the averaging principle (cf. chap. 10, §52] 
or jl2j ^ in order to simplify the effect of the perturbation on the three- 
vortex system. In future work we will deal with the delicate problem of 
determining the precise conditions under which the averaged dynamics is a 
good approximation for asymptotically large times as the binary becomes 
small. 

Let {zf, Ti}, i = 1, 2, 3, 4 be the complex coordinates and vorticities of a 
system of four point vortices such that initially the distance between vortices 
zi and Z2 is very small. The symplectic matrix in coordinates {xi,yi), with 
Zi = Xi + \yi, that corresponds to (jSJ, is given by 



Jo 
-Jo 



with J( 







V 



r4/ 



We start with a transformation Ti that takes the original complex coordi- 
nates to the coordinates {(", z, Z3, Z4}, where C, is the vector from zi to Z2 
and z denotes the position of the center of vorticity of the binary system 
(zi,Z2). Then Ti is the linear transformation given by the matrix TZ{Ti), 
where M : GL(n, C) — > GL(2n,M) is the operator transforming a complex 
matrix into its real form and 



/ -1 


1 









Ta 


















1 





V 








1/ 



A calculation shows that Jo Ti = Si with 

/ FiFa 

Ji 



-Ji 



and Ji 



+r2 



V 



TI + T2 



so that Ti is a generalized canonical transformation. 

The coordinates in the target space of T2 already correspond to the two 
systems in which we want to split the four-vortex problem. However, the 
dynamics in in the 3- vortex system needs to be reduced in order to simplify 
the description of the dynamics. To this end, following we introduce a 
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discrete Fourier transform given by: 



N 



Qn + ipn = ZaExp[i2im{a - \)/N] 



a=l 



applied to {z, Z3, z^} with = 3. The resulting linear transformation, tak- 
ing the canonical coordinates {x, xq, x^, X4,y, yo, 1/3, to the new canonical 
coordinates {x,qo,qi,q2,y,Po,Pi,P2} is given by T2 = U{f2) with 



T2 



\ 

1/^3 i/\/3 i/ys 

1/V3 e2'^/\/3 e-2i^/V3 

VO 1/^ e-2*^/V3 e^'^'/VsJ 



Let J2 be defined by T^Ji T2 = Jl2- A calculation shows that, J2 = 
A B 



^ I with A antisymmetric and B symmetric. In order to have a 

canonical transformation in the generalized sense we need to impose the 
condition ^ = and that -B be a diagonal matrix. A calculation shows that 
this happens if and only if 

Ti + r2 = Ta = r4 := r . (22) 

In what follows we assume this condition. A calculation shows that in this 

case J2 = Ji . 

A further reduction is achieved by introducing semi-polar coordinates. 
This is a standard procedure in celestial mechanics. The corresponding 
transformation taking the coordinates {x, qq, qi,q2,y,Po,pi,P2} to the new 
coordinates {j,jo,ji,j2,0,&o,0i,d2} is defined by 

F3{x,qo,qi,q2,y,Po,Pi,P2) = 



( 



x'^ + y'^ qo+Po Qi+Pi qI+pI 



2 ' 2 ' 2 ' 2 ' 
ArcTan(y/a;), ArcTan(po/go), ArcTan(pi/gi), ArcTan(p2/g2)) ■ 

A calculation shows that this is a canonical transformation with induced 
symplectic form J3 = J2- 

A final transformation is motivated by the fact that the first term Hq in 
the Taylor expansion of the hamiltonian in the variables {j, jj, 6, 9i}, taking j 
as the expansion parameter, depends on 61 and 62 only through its difference. 
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(Notice that Hq can be interpreted as the hamiltonian of the 3- vortex system 
unperturbed by the presence of the binary-system.) Thus we introduce the 
transformation T4 induced by the generating function 

G3 = -^j - ifoh + -j2)+ ^2{jl + ^2) , 

and 9i = —dG^/ji , ii = —dG^/dipi. The induced fourth symplectic matrix 
defined by T4 J3 T4 = J4 is given by 

\ 

r 

r 

V ry 

and with T as defined in ()22() . The transformation T4 takes the canon- 
ical coordinates {j, jo, ji, j2,G,6o,6i, 62} to the new canonical coordinates 
{i,io,ii,i2,ip,(fo,(fi,(p2}- 

5.2 Averaging and hamiltonian reduction 

We are now in a position to try to simplify the hamiltonian of the four- 
vortex problem using the canonical variables {i, io,ii,'i'2, f, Vo^ Vi; ^2}- The 
variable i represents the square of the radius of the binary-system. Let 
e = ^/2i and write 

H = Ho + eHi+ e^H2 + 0{€^) . 
A calculation shows that 

A^Hq = 2T2 {-T+T2) log(e)-r2 (log(-2 {i2+^-ii-i2Vn^2 cos(2(^i))) 
+log(\/-n -i2y/h- h cos(2(/3i)-2 {i2+V^^/-il - 12 - h cos((/?i) sin((/?i))) 
-flog(-2i2+V-n - 12 - i2 cos{2ipi)+V^yJ-ii - 12 - 12 sin(2(^i))) . 

Notice that the first term in the r.h.s. is proportional to log(e) and thus 
blows-up as e ^ 0. However, since Hq is cyclic in (p,ip2, we have that, in 
the unperturbed dynamics, e is, not surprisingly, a constant parameter. 

It turns out that the second term in the e-Taylor expansion of the hamil- 
tonian vanishes; that is to say, Hi = 0. The third term, however, is an 
important one and has a very complicated expression. To make it tractable 
we apply the averaging principle assuming that the ip, the variable describ- 
ing rotation in the binary-system, is a fast variable. Although we discussed 



J4 



J4 
-J4 



with J4 
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above conditions under which this is a reasonable assumption, a rigorous 
study of the time scales in which can be considered a fast variable is left 
for a forthcoming paper. 

The averaging principle is applied formally by space-averaging the hamil- 
tonian H2 with respect to the variable ip, i.e. 

1 Z"^'' 

H2 = — H2dip. 

2vr Jo 

A calculation implemented in Mathematica yields that 

fv fiiH,i2,^i) 

where 

/i := (-l + 27r) T2^ {-T + T2) {si2^-l2^-ii - is vV^i2^ cos(299i) 
+ 6z2 ( - ii^ + 12^) cos(4<y9i) + ii^ V-n - «2 - ^2 COs(6(/9i) 

- - 12 Vh - ^2^2^ COs(6<y9i)^ 

and 

/2 := 27rr ^ii^ + 3^2^ - 4 V-ii - ^2 - ^2^2 cos(2(^i) 

-2(ii2-i2^) cos(4<^i))^ 

Notice that H2 is again cyclic in (p2- (We remark that this is not true for 
the original H2.) Hence, from the averaging principle, we obtain a family 
of one-degree-of-freedom systems, parametrized by the constant parameters 
6,^2, having a two-dimensional phase space in the coordinates {^1,951} and 
governed by the hamiltonian H := i?o(*i) ¥'1; 12) + £^-^2(^1, y'l; ^2)- This is 
interpreted as describing the dynamics, in the averaged approximation, of a 
3- vortex problem where one of the vortices is a binary of radius e. 

6 Conclusions 

We have exploited the (degree minus one) homogeneous character of the 
equations of motion in the square distances of the four-point vortex problem, 

together with energy and inertia considerations, to obtain various character- 
izations of candidate evolutions leading to partial or total collision. Also, for 
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the general AT-point vortex problem, we showed that regular total collisions 
imply that the virial of the system is zero. Motivated by the relationship 
between candidate total and partial collapse we began the study of a four- 
vortex system containing two vortices close to binary collapse using the 
method of hamiltonian averaging. A detailed study of the conditions and 
time scales under which the averaged approximation remains valid is ma- 
terial of a forthcoming paper. In that work we plan to also generalize the 
reduction procedure described in this paper. 
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